function exp = taylorsimple(A)
exp=zeros(size(A));
add = eye(size(A));
for i = 1:1000
    exp = exp+add;
    add = add*A/i;
    if(norm(add)==0)
        break
    end
end